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We study the Z(4430)* meson within a heavy hadron molecule interpretation where the Z is 
considered as a bound state of a vector Z)*(2010) and an axial- vector Di(2420). We give predictions 
for the strong hidden-charm Z(4430)* -^ ir^ip and n^ip' decay widths and also study the radiative 
Z [J^ = 1~) —^ TT 7 decay properties in a phenomenological Lagrangian approach. Our findings 
are qualitatively in line with the experimental observation that the Tv^ip' transition dominates over 
the TV ip decay mode despite a smaller phase space. The width for the radiative mode tt 7 is sizable, 
^^ ' allowing possible detection in future experiments. 
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I. INTRODUCTION 

The spectrum, production and decay properties of observed hadrons still pose major challenges in their theoretical 
understanding. In particular, already in the light meson sector we have experimentally observed resonances, as for 
instance the low-lying scalar mesons, which cannot be simply and consistently explained by a leading order quark- 
antiquark (qq) structure. But an unique interpretation of such anomalous meson resonances is also not at hand, 
not even for a single case. Detailed discussions concerning the possible non-qq nature of meson resonances find their 
repetition in the heavy meson and especially the charmonium sector [l|, y| . Here investigations are essentially fueled by 
the enormous progress on the experimental side with the discovery of many new charmonium-like states by the Belle 
and BABAR Collaborations at the _B-factories (see e.g. recent conference proceedings 0]), but also in experiments 
by the CDF, DO and CLEO collaborations [j|. Properties of these newly observed X, Y and Z mesons cannot be 
easily explained within the standard charm-anticharm {cc) assignment. For instance, besides the overpopulation in 
the cc spectrum due to the numerous occurrence oi X ,Y and Z mesons decay patterns of these mesons arc in some 
cases in contradiction to the standard charmonium predictions. Here one example is the sizable ujcf) hidden charm 
decay width of the F(3940) [51-13 which should be significantly suppressed in the charmonium picture. Some of these 
new charmonium-like states are considered good candidates for possessing a hadronic substructure which goes beyond 
the standard cc assig nment ranging from quark-gluon hybrid mesons [Sl, l9| and tetraquark states [10| to dynamically 
generated states [111 or bound states of two mesons called hadronic molecules [3, ll^l ■ A review on the experimental 
situation with a first overview of the present theoretical understanding is e.g. given in [l|, 12, Isl, l6|. 

The observation of the charged Z(4430)^ by the Belle Collaboration [lj| presents so far the culmination of the 
series of newly discovered X, Y and Z mesons. The Z~^ was originally found as a relatively narrow enhancement in 
exclusive B -> KTr+ip' decays with mass mz = 4433 ± 4(stat)±2(syst) MeV and width Tz = 45t}3(stat)j:i3(syst) 
MeV. The BABAR Collaboration also searched for the Z+ resonance in B^ -^ Kn^ip' decays but without positive 
evidence for a narrow Z^ — > TT^ip' signal |l4{ . They report a branching- fraction upper-limit for example for the 
process B{B° -> Z(4430)"A'+, Z^ -^ V'^r") < 3.1 • 10~^ at the 95% c.l. It was speculated that the resonant 
structure, observed by Belle, might arise from interference effects in the Kt: rather than in the Trip' channel. After 
this report the Belle Collaboration |l5| reanalyzed their data sample in a full Dalitz-plot formalism now with a 6.4 
a signal for Z+ -^ ir+ij}' with mass mz = 4443ti2ll3 MeV and width F = IQltf^tf^ MeV. This rcanalysis confirms 
and supersedes the previous Z+ resonance parameters of [Ij]. The larger errors contain e.g. the uncertainties in 
the spin assignment of the Z^ (J = 0, 1) and in the orbital angular momentum in the B decay. The reanalysis by 
Belle results in a product branching fraction B{B^ -^ /v-Z(4430)+)S(Z(4430)+ -^ 7r+^') = {■i.2tllt\l) ■ 10"^ [II 
which is consistent with the BABAR upper limit. The decay mode ^/j'tt"*", assuming standard conservation laws, leads 
to an identification of the Z~^ as an isotriplet state with positive G-parity. The J^ quantum numbers remain to be 
determined. As in other cases of the X, Y and Z mesons, sec e.g. the 1^(3940), the Z+ also shows a sizable coupling 
to the hidden charm decay channel. As a consequence the partial decay width of Z+ — >■ 7r+'0' is expected to be on the 
MeV scale [l| . In comparison, for the conventional cc configuration open charm decay modes arc dominant whereas 
hidden charm decay channels are Okubo, Zweig and lizuka (OZI)-suppressed and typically result in decay widths of a 
few keV only [a, l6|. Since this resonant structure was observed in the invariant mass of tt -0' the Z as a charmonium 
state with isospin / = 1 is a truly exotic resonance which might offer the possibility to uniquely pin down a non-cc 
structure. A further key issue concerns the observation mode of the Z+. While the Z+ was seen by Belle in the "0' 
channel no signal was reported in the 7r+0 mode implying a large ratio R = F(Z'+ — >■ tt~^iJj')/T(Z'^ — > 7r"*"0). Note 
that the suppressed mode 7r+0 is favored by phase space. Although some arguments [l6l - ll9j were put forward to 
understand this dynamical selection rule, a full quantitative explanation is not given yet. The Belle Collaboration 
reported [20j on further charged states in the charmonium sector. There is evidence for two charged resonances in 
the TT+Xci channels termed Zj^(4050) and .^2~(4250). However, the signal is much poorer than for the Z(4430)+. 

Several non-cc structure interpretations have already been discussed in relation to the Z(4430)^ (since the Z^ carries 
charge a cc hybrid configuration is obviously excluded) . The Z^ is considered a candidate for a radially excited ccud 
tetraquark [ifl: ES US llll , for a less compact hadronic meson molecule [22,123 o^ be just due to threshold effects [19| . 
Further explanations for the occurrence of the Z(4430)+ are a cusp in the DiD* channel [2J|, a radial excitation of a 
cs configuration j25| . a baryon-antibaryon (baryonium) bound state [26[ or even a ip' bound state in mcsonic matter 
[27| . The molecular interpretation is rather natural since the Z(4430)"'" mass lies extremely close to the _Di (2420)1?* 
threshold at 4.43 GeV. Assuming an S'-wave Di{2A20)D* bound state possible quantum numbers are J^ = (0, 1, 2)~. 
This would also imply that the Z^ — > ir^ip' decay proceeds in a final state P-wave as opposed to an S'-wave for the 
J^ = 1+ tetraquark proposition. Note that the I?i(2420) is fairly narrow with a total width of about 20 MeV, which 
would principally allow the formation of a hadronic bound state. Alternatively, the bound axial- vector charm meson 
can also be identified with the D[ = _Di(2430) which within errors is degenerate in mass with the I?i(2420) but has 
a rather large width of about 400 MeV. 

Aspects of the hadronic-moleculc interpretation have been studied within various models including meson-exchange 



potential approaches [28l - l30| . QCD sum rules [SlllS^I and cfFcctivc Lagrangian techniques [18|, |33[. First analyses [28|, 
^^ based on the long-range one-pion exchange (OPE) mechanism conclude that the isovector DiD* or D'^D* systems 
cannot form a J^*^ = 0~~ or 1~~ bound state. Based on the binding energy of the DiD* system of around a few 



MeV the OPE study in [28| suggests a probable isovector 1 ^ assignment for the Z+. Further inclusion of sigma 
meson exchange leads to S'-wave binding for D[D* with J^ = 0^, 1^, 2~ [SO], but the large width of the D[ probably 
disfavors the formation of a molecular state. Inclusion of the sigma-exchange potential also leads to binding for the 
DiD* configuration but only for J^ ~ 0~ and at the price of a large cutoff which in turn leads to an enhanced 
attraction of one-pion-exchange. An evaluation in the context of QCD sum rules [3^ favors the molecular DiD* 
bound state interpretation of the Z+ with quantum numbers 0~. The study of the low-energy DiD* interaction in a 
quenched lattice calculation also indicates attraction in the J^ = 0~ channel, but this effect is considered probably 
to weak to lead to the formation of a bound state [Sa] . 

Two-body decays of the Z+ were also analyzed in effective Lagrangian methods. The open charm decays D+D*", 
]j*+]jo and D*'^D*'^ were analyzed in [3j] and, in spite of ill determined coupling constants, argued to be suppressed 
in the molecular interpretation while dominant for the tetraquark configuration. More importantly, the Z^ was 
observed in the hidden charm mode tt"*"-;/'' while no signal was seen in the tt "*"')/; decay channel. In the unpublished 
work of [l^l these channels were investigated in addition to the dominant D*D*tt decay, but the conclusion on the 
possible suppression of J/tjj depends very much on form factors and the regularization in the loop diagrams. 

In the present work we reconsider and pursue a quantitative explanation of the hidden charm decay modes Z^ — > 
TT^'ip and 7r+-0' in the context of a molecular DiD* bound state interpretation. In addition we determine the radiative 
decay width Z+ -^ 7r+7 as a further key feature of the molecular idea. As suggested by the above mentioned studies 
related to possible binding we consider the quantum numbers J^ = 0~ and 1~ for the hadronic molecule. In technical 
aspects we proceed in an alog y to the open and hidden-charm hadrons £'*o(2317), Dsi{2460), Ar(3872), y(3940), 
y(4140), A(2940), etc. [tI. [36l - l39l | considered previously as hadronic molecules. For our analysis we use an efi'ective 
Lagrangian approach for the treatment of composite objects — molecular states [TI. I36l439l | . The hadronic bound state 
is set up by means of the compositcness condition |40| which also allows for a self-consistent determination of the 
coupling strength between the hadronic molecule and its constituents. 

This work is organized as follows: In the following section we discuss the set up of the mesonic bound state and 
introduce the effective Lagrangian approach which we use to study the decay properties of hadronic bound states. In 
the subsequent sections we apply our method in order to compute the radiative Z+ — > 7r+7 decay in section lllll and 
strong Z+ — >■ Tr~^ ip^'^ decays in section ITVl Our results are presented in section [V] At the end of this work (in section 
IVip we give a short summary of our findings and draw the conclusions. 



II. THEORETICAL APPROACH 

In the present study we assume the Z~^ to be a pure bound state of an axial D^' and a vector D* meson. In the 
charmed meson spectrum two nearby P-wave excitations with J^ = 1+ are expected. These two axial Di states 
can be identified with the Z?i(2420) = Di and the Di(2430) = D'l. In the heavy quark limit the two degenerate 1"*" 
states are characterized by the angular momentum jg of the light quark with jq = 3/2 and 1/2. While the strong 
decay Di{jq = 3/2) -^ D*n proceeds by D-wave, the transition Di{jq = 1/2) — >■ D*tt has a final S'-wave. The state 
decaying via D-wave is narrow while the one decaying in an S'-wave is expected to be broad. Since heavy-light mesons 
are not charge conjugation eigenstates the axial states can also be written as a superposition of the ^Pi and ^Pi 
configurations {J = L and S = or 1) with 

\Di)= cos0|Vi)+sin0|Vi), 
\D[) = -sin(/)|Vi)+cos0|Vi). 

More detailed analyses |4l| - t43j of the mixing scheme in terms of the total width indicate that the mixing angle has 
a value of about (f> — arctan (l/v^) ~ 35.3°, the "magic" value expected from the heavy quark limit. With this 
phase convention (alternatively (p = — arctan(v2) ~ —54.7° can be used) the Di state is identified with the narrow 
1^1(2420) (r sa 20 MeV (H), while the broad i:)i(2430) (F w 380 MeV jH) is connected to the D[. Since our aim is 
to study the Z+ as a mesonic bound state containing a Di state the narrow _Di(2420) with its long lifetime is more 
favorable than the broad Di(2430). 

Since the ^(4430)+ was observed in the ■0'7r+ final state, isospin and G— parity assignments are I^ = 1+. If the 
2'(4430)"'" is a S-wave DiD* molecule the J^ quantum numbers are 0~, 1~ or 2~. Here we restrict the study to the 
0~ and 1~ cases since J^ — 2~ seems excluded by the small phase space in the B — >■ Z~^K production process |45| . 



Following the convention discussed in [29|, |45| the particle content of the isospin multiplet is given as: 

1 



\Z+) = —{\DtD*^) + \D\D*+)) , 

1^0) = \{\^B+D*-) - \d'(D^) + \D~D* + ) ~ \l5lD*^)) , (2) 

\Z-) = -^^{\D',D*-)+\D^D*")). 

In the present method the meson bound state is first set up by the effective interaction Lagrangian between the 
hadronic molecule and its constituent mesons. In case of J^ = 0^ the Lagrangian reads 



C 



ZD,D' 



z-(..)/dy$(j/^){i?r(:^- !)<(-+ !)+<(-- f)<(^+f)}+h.c.. (3) 



V2 
For J^ = 1^ the respective Lagrangian Lzd^D' is given by 



with X and y being the center of mass and relative coordinates. The correlation function $(y^) takes into account 
the distribution of the constituent mesons in the Z resonance. Its Fourier transform enters as a form factor in our 
analysis which also leads to a regularization of the structure integrals. In earlier works [4y, |43] it was found that 
obscrvables as decay widths which are studied in this paper are not sensitive to the specific shape of this form factor 
as long as the intrinsic scale remains the same. In the present case we deal with the Gaussian function 

'^^y'^ = l-0^^''''^^~'''^^ $(fc|) = exp(-fc|/A2), (5) 

where the index E refers to the Euclidean momentum. The size parameter A is fixed in the physically meaningful 
region of a few GeV. In the present work we study finite size effects by varying Az in the range of 1.5 — 2.5 GeV. In 
the special case of a pointlike molecular structure the correlation function $(y^) reduces to the delta function 5^{y), 
equivalent to <&(fc|;) — >■ 1. 

The coupling of Z to the virtual constituents, denoted by gzDiD* , is fixed by means of the compositeness condition [7|, 
l36l - [40l . l4g| which provides a self-consistent method to adjust this quantity. Since we deal with a hadronic bound state, 
the physical field of the Z meson should be fully expressed by the fields of the constituent mesons and therefore does 
not contain a bare component. This is achieved by setting the field renormalization constant Zz to zero with 

Zz = l-9l,^„,t'{ml)^0, (6) 



where n'(m|) = g^ ^'E,'{mz) is the derivative of the mass operator illustrated by the diagram in Fig. [T] If the 

Z meson is a vector-like object (quantum numbers J^ = 1~) the mass operator is split into its transverse and 
longitudinal parts S and E^ 



E^^ip') = I](/).gr + S^b^)^ , (7) 

where g'^ = g^"^ — ^-^ and g^p^i = 0. In the compositeness condition of Eq. © only the transverse part enters. 




FIG. 1. Mass operator of the Z(4430) meson 



Meson loop diagrams are evaluated by using the free meson propagators, which in momentum space read as 

in case of pseudoscalar and scalar mesons {Hi = P,S). For the case of vector and axial- vector mesons {H2 = V,A) 
we use 

Here we do not include the finite widths of the constituent mesons in the propagators (especially for the -Di(2420) 
with a total width of about 20 MeV) since this effect is negligible in the present analysis. 

The mass values for the intermediate and final state mesons are taken from [4J| . For convenience we also introduce 
the binding energy e defined by the difference between the central value of the Z+ mass and the lower threshold 
{D^W") with 

mz+ = m^+ + m-^, — e. (10) 

We finally present our results in dependence on the possible values of the binding energy e. 

III. RADIATIVE DECAY 

We first consider the radiative decay Z+ — > 7r+7 which in the molecular DiD* interpretation proceeds by the 
diagrams of Fig. [2] In the following we discuss the relevant vertices entering in the radiative decay process. 

One of the interaction Lagrangians relevant for the radiative decay arises from gauging the free Lagrangians by 
using minimal substitution 

d^M^ ^{d''TieA'')M^ , (11) 

which leads to 

£^"(1) = ieA^{Z-d^'Z+ + 7r-a%+ + g°'^V-d'V+ + 9^'^V^d"V+ - g^'°'d'^V-V^) , (12) 

where Ad^B ~ Ad^B — Bd'^A and V ~ D* , Di. The corresponding vertices enter in Figs. [11(a) and[2](b), where the 
triangle diagrams are obtained by coupling the final states to the Z'^ constituents. These graphs yield the dominant 
contributions to the decay amplitude. Because of their nonlocal structure the strong interaction Lagrangians ([3]) 
and ([4]) are not invariant under C/em(l) transformations and need to be modified accordingly. We use the method 
suggested in [49| where each charged meson field M is multiplied by an exponential containing the gauge field 

M^{y) -> eT'^(^'"'^)M±(j/) (13) 

V 

and I{y,x,P) = J dZf^A'^{z). This modification leads to further interaction vertices contained in Figs. [2] (c) and 

X 

[2] (d) . These additional graphs [2] (c) and [2] (d) are strongly suppressed but they have to be included to guarantee full 
gauge invariancc. 

The interaction between the final pion and the charmed mesons Di, D* in the loops is set up by the interaction 
Lagrangian 

CD,D'.^^^^D'^''7TTW^, + h.C., (14) 

where V^"^ = d^^V^ — d'^V^ is the stress tensor of the vector mesons V = "07^1 ^^^ D* ■ The interaction vertex 
involving the I?i(2420) meson should contain a dominant D-wave D*7r coupling (see e.g. [4l|, |42|, [SO] ) which dictates 
the form of the Lagrangian in (fT4|) with two derivatives involved. This dynamical selection rule, obtained in the heavy 
quark limit, also leads to the form of the Lagrangian ([22]) in case of the strong decays discussed in the next section. 
The coupling constant goiD'-n is derived from the width of Di — >■ Z?*7r which is the dominant decay mode of Di. The 
partial decay width is expected to be around 20 MeV Hil, where T{D\ -^ D*+n-) = 2r{D"i -^ D*°tt°). The decay 
width is set up as 



IGttto? 



T{Dl^D*+Tr-)^ ' ^1' f "' 1X1, (15) 



where A(a, 6, c) = a^ + 6^ + c^ — 2ab — 2ac — 26c denotes the Kallen function and \M.\ = 1/3 X) \-^\ represents 

pol 

the spin-averaged and summed over polarizations transition amphtude squared. The effective i'-wave interaction 
Lagrangian of Eq. ([T4l) leads to the matrix element 






Hence, the decay width for a D-wa.ve decay is of the form 



9DiD-7,[PD-PDig^ -Pd'PdJ 



T{D^^ -^ D*+7r 



-N _ Qp^D'Ti .1/2 



OGttttt, 



l^AV2(A + 6ml,^m|,.) 



with the resulting coupling gD^D*■^ = 0.49 GcV^^ for T{Di -^ D*tt) «20 MeV 

Q 






(16) 



(17) 



(c) (d) 

FIG. 2. Diagrams contributing to the radiative Z'^ -^ 7r"'"7 decay. 

In case of J^ = 0^ the radiative decay Z+ — > 7r+7 is excluded due to gauge invariance. If the spin-parity of the 
Z+ is 1^ then the transition amplitude Ai'^'^ has the structure 



M''''^eFz^^iml,mlO)e''^^''p, 



qp: 



(18) 



where gzTvy is the effective coupling. It is related to the corresponding transition form factor gzn-y — FzTT-y('ni%,m'^, 0) 
evaluated via the loops of Fig. [51 Finally, in terms of the effective coupling gz-n^ the decay width r(Z+ -^ 7r+7) is 
given by 



a 
24' 



nZ^ ^ -+7) = i,gl.,^% 1 - ^)' - ^,9l.,m% 






(19) 



IV. STRONG HIDDEN-CHARM DECAYS 



In order to study the strong hidden-charm decays Z+ — )■ t/i^'V^ wc first set up the interaction between the final 
state and the constituent Di and D* mesons. Wc use the effective Lagrangians 



Cd'Dit 



9d* Dtt 7-)* t 



D*,^TTTid''D + h.c.. 



2V2 "^ 



c 



DiDip 



g D^D 



"/' n'^'^ 



D^'^iP^.^D + h.c. 



(20) 
(21) 
(22) 



The coupling of the Z+ to its constituents is determined by Eq. ^. In addition, CdiD'-k was already defined in 
p4[) in the framework of radiative decays. The coupling strengths gu'D-n — 17.9 and gD*D*j/ip ~ 8 as well as the 



ratio gtp'D'D' /g-tpO'D' arc taken from heavy hadron chiral perturbation theory (HHChPT) [3a, |3g. This ratio of 
couphngs of excited i/j' to ip is fixed by g^'O'D" /gipD'D* = tri^i f^/(Tn^f^i) = 1.67, where f^u) is the leptonic decay 
constant. An estimate for the remaining couphng gp, n ,h is taken from a coupled channel analysis (see e.g. [5l[). We 
use |5DiD*7r| ~ 72 MeV and \gDiDip\ ~ 29 MeV from [531 i^i order to estimate the ratio of the coupling strengths 
I"! = goiDipl goiD'-K ~ 0.4 ± 0.2, where we assumed an uncertainty of 50%. By using the ratio ri and goiD-K = 0.49 
GeV~^ as defined above, we can give a rough estimate of goiDit ~ 0.2 ± 0.1 GeV"^. As discussed below, the diagram 
(a) of Fig. O clearly dominates the hidden-charm decay. Therefore, in leading order the .Z"*" — > ■)/)(') tt"'" transition can 
be regarded to be proportional to guiD-4i^')- ^ variation of .to^D^w modifies the decay width accordingly. In order 
to estimate the ratio of couplings r2 = goiD^i' /goiD^i we use the "^Pq niodcl [53j . where the details of this procedure 
can be found in Appendix [^ Wc find for the ratio of couplings gDiDtp' I gDiDtp a value which is close to 2 and of the 
same order as the above mentioned ratio g^p'O'D' / g-ipO'D* = 1.67. This hierarchy of couplings involving the ^ and '0' 
charmonium states is consistent with the HHChPT scaling indicated above. Further on we include uncertainties in the 
predictions of the ^Pq model which for example might arise due to variations of the quark pair production amplitude 
(as e.g. discussed in l54i ) which is usually fitted to data. In literature the quark-pair production strength ranges 
between 0.4 [41[ to 0.5 [53. We therefore consider an uncertainty of 50% in the ratio r2 = goiD-ii' / goiDii = 2 ± 1. 
The diagrams describing the hidden charm decay are illustrated in Fig. [3] 





(b) 



FIG. 3. Diagrams contributing to tlie Z — >■ V' tt decay. 



Provided that the Z^ is a pseudoscalar the transition amplitude for the hidden charm decay mode can be expressed 
by two form factors Pi, 2: 



M%._ 



,^+^ = Fi{m%, ml,ml)p'^ + P2(m|, ml,ml)p'^ . 
Here only the first form factor contributes to the decay width T{Z^ -^ '0'')7r+): 



r(z+ ^ v''^7r+) 



•^^''^ 77712 2 >'Hml,ml,ml) 

LM =g 5-^-^5 

167rm| ' ' '''"* 647rm|m^ 



(23) 



(24) 



where g^^j, = Fii^Z' "^^i '^^) ^^'^ M^^ V^ z) ~ x^ + y^ + z^ — 2xy — 2xz — 2yz is the Kallen function. 
If we deal with a vector Z^ the matrix element is given by 

M"" = P3(m|,m|,m^)e"^'^VaP/3 
and by analogy we use g^^^, = F^{rn\,rn\,rn'^) in order to calculate the decay width 

,2 ^h'm%,rnl„ml) 



r(z+ ^ Vf''7r+) 



967rTO| 



(25) 



(26) 



RESULTS 



In the following section we discuss our results on the decay properties of the Z^ . As far as the hadronic molecule is 
involved the coupling constant is fixed by the compositeness condition. The corresponding couplings are summarized 
in Table H] for the quantum numbers J^ = 0~ and 1~. Values for the coupling constants are presented in dependence 
on the binding energy e. Since the Z+ mass lies close to the (DfD*^) threshold at 4.43 GeV the values for e are of 
a few MeV, in particular we vary e between 1 and 10 McV. The errors on the numerical results are due to variations 
of the model parameter Az from 1.5 to 2.5 GeV. 



TABLE I. Coupling constants gzDiD- (J^ = 0") and g'zDiD- {J^ = 1") in GeV for Az = 1.5 - 2.5 GeV and e=l - 10 MeV. 



[MeV] 



10 



gZDiD' {J'' = 



3.8 ±0.1 



5.6 ±0.1 



6.8 ±0.2 



g'zD.D' {J'' = r 



1.2 ±0.1 



1.8 ±0.1 



2.1 ±0.1 



TABLE IL Decay widths rz+^^+ in keV for J^ = 1" with e = 1 - 10 MeV and Az = 1.5 - 2.5 GeV. 



e [MeV] 



10 



■ Z+^ir + ' 



[keV] 



0.3 



+0.2 



0.6: 



0.8: 



The radiative decay width for Z+ -^ 7r+7 is analyzed for 1^, whereas for J^ = 0^ it is forbidden. The resuhs 
are indicated in Tab. |lT] for different values of the binding energy and the size parameter. Our predictions are in 
general rather sizable and of the order of 0.2 to 1.2 keV. In Tab. |ll]the smaller value of each entry corresponds to 
Az = 1.5 GeV while the larger one is related to Az = 2.5 GeV. 

The decay widths for the hidden charm decay channels Z — >■ Trip and Z -^ nip' are given in Tab. Illll for J^ = 0~ and 
in Tab. IIVI in case of 1~. In both tables the error bars indicate changes in the finite size effects with Az varied from 



TABLE III. Decay widths T^ 



iilCI; 



in MeV for J^ = 0", Az = 1.5 - 2.5 GeV and ra 



■ gOiD-ip' IdDiDi, 





e [MeV] 


Tz^^i, [MeV] 


r 

r2 = l 


z^^^> [MeV] 

ra = 2 ra = 3 


7? = r 

r2 = l 


Z^TTtp' f^Z^TTil 

ra = 2 ra = 3 




= 0.2 


1 


o.2+;!:i 


1+"° 


■^+°i 


0.81^:^3 


«0.2 


« 1.5 


f»3.8 




5 


0.41-j 


1+"-^ 


0.6 ±0.2 


i-5t^:^ 


«0.3 


f5il.5 


f«3.9 




10 


0.5+°-J 


O.lt^o 


0.81?,:^ 


n Q+0.5 


f»0.3 


^ 1.5 


f»3.7 




= 0.4 


1 


0.9tll 


o.3l^;l 


i-sir* 


q r+0.9 


fa 0.3 


«1.6 


fa 4.0 




5 


1 7+0.5 
J^-'-0.6 


0.61-J 


O n + 0.6 
^•0-0.8 


6.7tli 


f«0.3 


« 1.7 


f»4.0 




10 


n q+0.7 
^"J-O.S 


o.8l°:^3 


3.7i?:? 


o 7+1.8 


«0.3 


«1.6 


f«3.8 




= 0.6 


1 


2.11°;? 


o.8l[!J 


q t;+0.9 

'J-'J-o.g 


8-3_2'.o 


RiO.4 


« 1.7 


f=i4.0 




5 


3.9lii 


1 t^+o-* 

J--O-0.5 


6.7li:? 


15.6+4:1 


fa 0.4 


«1.7 


fa 4.0 
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5.3li:? 


2 0+"-^ 


« 7+1-8 


20.ltli 


f=i0.4 


« 1.6 


f=i3.8 



TABLE IV. Decay widths F 



z^iiC. 



in MeV for J^ 



1 , Az = 1.5 - 2.5 GeV and ra = goiDii' /goiDi, = 1-3. 





e [MeV] 


Vz^^i, [MeV] 


Lz 
ra = 1 


ra = 2 ra = 3 


i? = r 

ra = 1 


Z-»7rt/.'/r 

ra = 2 


Z— >-7ri/' 

ra = 3 




= 0.2 


1 


0.1 


0.1 


0.31^1 


0.81°:? 


^QA 


«5.5 


^ 14.0 


QDiDtp 


5 


0.1 


0.1 


o.6±n 


1 ^+"-2 


^Q.^ 


«5.3 


^ 13.6 




10 


0.1 


0.1 


7+°'^ 

U-'-O.l 


1 q+0-3 

J^-J-O.4 


«0.8 


«5.2 


f« 13.6 




= 0.4 


1 


o.4l°:l 


0.3-0.4 


1 6+"-^ 

J^-D-0.4 


3-8_o:8 


«0.7 


«3.5 


f«8.3 


aDiDi, 


5 


0.8l°o-2 


o.6l°o:l 


9 0+0.4 
^■^-0.7 


6.911:^ 


^O.S 


«3.6 


«8.6 
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l.l+?,"3 


o.7tS:i 


o c+o.e 

■J-O-0.7 


8.rtl:^ 


«0.6 


«3.3 


f^7.9 




= 0.6 


1 


l 9+0.3 


n s+O-2 
0-8_o.i 


q Q+0.7 

•J-O-0.8 


9.01^:1 


«0.7 


«3.1 


f»7.4 




5 


2 3+°-* 


1 ^+0-2 

J^-O-0.3 


6.9±l;°e 


16.2i-« 


f«0.7 


«3.0 


f»7.0 
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9 q+0.6 


1 q+0-3 


8.7tlj 


2QA+Jd 


«0.7 


«3.0 


f«7.0 



1.5 to 2.5 GcV. The ratios ri ~ snioi, _ q^^ _|_ q 2 ^^^ „ _ u-^ny = 2 ± 1 wc fixed from values obtained in eoupled 
cliannel analyses 51| and from the phenomenology of the '^Pq niodel (see Appendix). We consider uneertainties in 



the ratios ri and r2 and study the decay properties of the Z^ for possible values of the coupling constants. 

The hidden charm decay mode is generated by the graphs of Fig. |31 where diagram |3] (b) clearly dominates the 
transition amplitude by one order of magnitude in comparison to grapli[3] (a). This dominance results from the sizable 
£)*Z?7r-coupling but also from the lighter D-meson mass (compared to D* exchange) in the rescattering process. As a 
consequence the decay widths are very sensitive to variations of the couplings g^j-^^i^ii) ; variations of ri and r2 enter 
approximately quadratically in the decay widths. For instance, Tz->-it,Ii' f^z-^Trtii oa r| as can be read off froin the 
results in the last three columns of Tabs. IIIII and ITVl 

For values of the ratios in the region ri = gDiDtpl goiD'-K ~ 0.4 to 0.6 and r2 = goiDip' / gDiDih = 2 — 3 the decay 
width Tz^-K-4)' is in the MeV range which is consistent with the expectation from observation [l|. However, if the 
ratio ri is relatively small as in the case of ri = 0.2 the decay width Tz^Trtjji becomes smaller than one MeV, which 
seems excluded by observation. 

For the quantum numbers J^ — 0^ the decay mode Z — ^ tt-^ is suppressed relative to Tr-f/;' by a factor « 2 for 
r2 = 2 and about 4 for ^2 = 3, respectively. In case of 1~ the ratio of decay rates R = r(Z+ — > 7r+-0')/r(Z+ — > tt+^/i) 
is even larger with i? « 3 for r2 = 2 and _R w 8 for r2 = 3 which is at least qualitatively in line with the experimental 
observation that the 7r+?/;' decay mode dominates the tt+'j/' partial decay width. The ratio R is rather insensitive to 
variations of the model parameter Kz and the binding energy e as indicated in Tabs. IIIII and IIVI for J^ = 0^ and 
J^ = 1~. On the contrary, for equal couplings goiDii' and goiD^p, i-C ^2 = 1, the branching ratio R = r{Z^ — > 
T:^ip')/T{Z~^ — >■ TT "'"'(/') is smaller than one which presently is in contradiction to experimental observations. 

At least in the molecular scenario, the radiative decay might help to partially settle the J^ quantum numbers of 
the Z, the strong Z+ — )• 7r"*"-0' decay is not very sensitive to the choice of J^ . The dependence on J^ of the strong 
hidden-charm decay widths is not very pronounced, for J^ = 1~ decay widths are only slightly smaller than for 0^. 

Within the hadronic bound state interpretation hidden-charm decays were also discussed in [1^ . But in the present 
evaluation the contributions of the D* and D rescattering processes represented in Fig. |3]are significantly different. 
In our case diagram [3] (b) is dominant because of the sizable coupling of the D* meson to the Dtt mode. On the 
contrary in [l^ the diagram in Fig. [3] (b) is highly suppressed by a factor of about 15 — 30 compared to the D* 
exchange process, which was explained by the small gDiDip coupling. We want to remind that the prediction for 
R = r(Z+ — > 7r+f/;')/r(Z+ — > 7r+^) primarily depends on the explicit values of the coupling gDiDMil ^'^d therefore 
on the ratio goiDtp' / gDiDijn while variations of Kz and e only play a minor role. In comparison, in |18| sizable values 
for the ratio R are only obtained for small hidden charm decay widths. 



VI. CONCLUSIONS 

In the present work we consider hidden charm and radiative decays of the Z(4430)+ in a Di(2420)£)*4-/i.c. molecular 
structure interpretation. As guided by previous studies of possible binding mechanisms in this system we choose the 
preferred J^ = 0^ and 1^ quantum numbers for the Z^ . 

In the predictions for the decay widths we study the influence of finite size effects and the dependence on the exact 
value of the binding energy. We give a first prediction for the radiative decay width ^(4430)+ — >■ 7r+7 only allowed for 
J^ = 1^, which is about 0.5 to 1 keV. This order keV result should allow for a possible detection. We also analyzed the 
hidden charm decays Z+ — > tt"''?/', presently not observed, and Z+ — > ir'^ip' , the discovery mode of the Z'^ . Both decays 
are generated by D* or D rescattering processes. Predictions depend crucially on explicit values for the couplings 
9DiDib^''' ■ For approximate values of these couplings we took guidance from a coupled channel analysis |5ll.l52| (which 
essentially fixes goiDi)) and in addition from the ^Pq model (for the ratio r2 ~ goiD^i' /goiDip)- Low values for r2 ~ 1 
generate predictions for Tz^ttiP' which are much to small to justify observation and for R = Tz^nti)' /^z-^ttiP which 
contradicts the nonobservation of the tti/j mode. Intermediate and larger values for r2 « 2 — 3 lead to order MeV 
predictions for r^-^Tri/'' and also to a ratio i? sa 2 — 9 at least qualitatively in line with experimental constraints. For 
J^ = 1~ the ratio R is slightly larger, but in general a sizable dependence of the predictions on the choice of J^ 
is not observed. Also, finite size effects (as contained in the size parameter Az) and the exact value of the binding 
energy do not have a large influence on the predictions. 

Present predictions have a slight tendency to support the DiD* bound state interpretation of the Z(4430)+ at least 
what concerns the hidden charm decay modes. A further evaluation of open charm decay modes does not seem to be 
decisive since predictions also depend crucially on principally unknown coupling constants j33| . In this respect the 
radiative mode 7r+7 offers a further test for the structure issue of the Z+. 

From the experimental side it is clear that the existence and, if possible, quantum numbers of the Z(4430)+ still 
have to be flrmly established. As long as the charged charmonium-likc structures are not ruled out they remain an 
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interesting object for future experiments e.g. the Z^ could be studied in nucleon-antinucleon scattering processes in 
the upcoming PANDA experiment j56| . 
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Appendix A: Ratio of couplings gDiDi/.'/gOiD'/' 



Here we give an estimate for the ratio of couplings for DiDJ/ip and DiD^p' by using the ^Po rnodel j4lLl53l.l55l.l57|. 
The ^Pq model is a standard phenomcnological tool to analyze hadron decays. Thereby a qq pair is created from 
vacuum with quantum numbers I'^{J^'-^) = 0+(0"'""'"), hence ^Pq in spin-orbit coupling. The ^Pq model is rather 
sensitive to variations of the parameters i.e. the elementary pair creation strength and the radii of the hadron wave 
functions involved. The model can deliver meaningful results for strong hadronic decay rates when evaluated in the 
center of mass frame and provided that all initial or final state particles are on-shell. Further extensions of the ^Pq 
model concern for exam ple nucleon-antinucleon annihilation processes |58l . |59[ and the determination of baryon meson 
coupling constants |4ll . |60| where the emitted meson is not necessarily on-shell anymore. Since in the present work 
we deal with transitions between off-shell mesons, the Di — >■ Di}} and Di — > Z?-0' decays are kinematically forbidden, 
we cannot fix the couplings directly from the model. However, since the ratio of transition matrix elements has a less 
pronounced parameter dependence, we use the ^Pq model to determine the ratio of the couplings r2 — goiDti'' / QDiDiji- 

In the ^Pq model the transition amplitude for the process A — > BC is given by 

T^_BC = (*f,,,,™,(l,3)vI'^^,^,™^(2,4)|03p„(3,4)K^,^^,„Jl,2)), (Al) 

where the indices i = 1, 2, 3, 4 refer to the respective quarks. For the initial and final mesons we use, as usual, simple 
harmonic oscillator wave functions 

^n,iAP) = A^„.i4+'^'(^V)exp( - ?^p^){Rp)i[\{siS2)sms)®Yi,,{p)\j^^^ (A2) 

x5^^\P-p^-p2)xln{^2)x^'^''{<li<l2), 



with the normalization iV,, / = ('_j'\2Tt+; / 2nR ^^^q radius R and the relative momentum p = mip2 m2Pi ^ 

Quark i is characterized by its mass rrti and spin s,;. Lj^ and Yim represent the Legendrc and the spherical harmonics, 
respectively, x'' ^nd X'^^'^ denote the spin and flavor-color wave functions. 

The additional quark-antiquark pair required for the decay of a meson into a two-body final state is generated by 
the ^Pq operator 

O.p, = XVs%U'^^^ (pi - pO<5(3) (p., - p,') (A3) 

spectator quarks 

where A is the constant qq production amplitude and 

V'po^ = E(-)'^''(l 1 - *" *" |00)X,(P3 - P4)S^'^ {Pz - P^W'f (A4) 

with the Pauli matrix a and yim{p) = |p|^im(p). 

In the following we give the final result for the amplitudes characterizing the ip^'^ — > DiD transitions (2)l'^5i — t-^ 
Pii^Pi) +^ S'o, where only the _D-wave contribution [L = 2) is relevant 

Tl=2 = \8^^\Pa -Pb- Pc) J2 YLmAP)Tr''yLil^,mi^,mL) 



with 



L,mL 



rpspacc,n=0 .VJlZ /I /, 2 n 2 n)2 / ^ n2 2\ 

i?2i [T/ll 7 o 2 . . , 4 . . o\„.. o , i?2 



and 



rTispace,n=l r,-^ I 2 ^ 3 d2 7-)2 , ^ 4 d2 n)2 \ d2 n)2 / ^' o2 2\ 

Tl12 =-3*^j-^^-(-a--a ^-aRP + -a R P )R P cxp ( - — P a ) 



h{fi,mi^,mL) ^ (110 0\2 0) (11 - A^ - m,^ |LmL)(-)™^Ji,2 (A5) 



(see also in e.g. [54|). Here L is the relative angular momentum between the final mesons and Jbc is the total spin 
of the two final states B and C. The different quark masses of the charm mesons are accounted for by the factor 



rrig+TTic 
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The latter expression (jASp represents the spin part of the amphtudc for a ^Pi and ^Pi state, respectively 



rfp" = y -—--{JBc2'mBCiTiL\lmA)SL.,2 , 
r|p" = a/-— ^(jBc2mBC™L|l"M)(5L,2- 



Finally we consider the mixing \Di) = \/||^^i) + yhl'^^i)^ which leads to 

Due to this mixing pattern the Di couples via £)-wave only. 

We find that the ratio goiDii' / gDiDip = T^'^fj^'"" jj^'^^^c.n- _ 2jg jg independent on the radii R and coupling 
strength A, where we evaluated the amplitudes at threshold i.e. P=0 which is in analogy to the determination of the 
nucleon-meson couplings in [5J]. The quark masses m, = 0.33 GeV and TOc = 1.6 GeV are taken from j4lj . 



